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1.  INTRODUCTION 

This  is  SARA’s  9h  Quarterly  Report  for  “Breakthroughs  in  Low-profile  Leaky-Wave  HPM  Antennas,” 
a  37-month  Basic  Research  effort  sponsored  by  the  US  Office  of  Naval  Research  (ONR).  This  work 
includes  fundamental  theoretical  analyses,  numerical  modeling,  and  related  basic  research.  Objectives 
include  to  discover,  identify,  investigate,  characterize,  quantify,  and  document  the  performance,  behavior, 
and  design  of  innovative  High  Power  Microwave  (HPM,  GW-class)  antennas  of  the  forward-traveling , 
fast-wave ,  leaky-wave  class. 


1.1.  Overview  of  Previous  Activities  (1st  thru  8th  Quarter) 

During  the  first  quarter,  we  prepared  and  established  useful  equations  and  algorithms  for  predicting 
reflections  and  transmission  of  incident  TE  waves  from  parallel-wire  grills,  dielectric  windows,  and 
combinations  of  wire  grills  with  dielectric  windows,  in  problems  reducible  to  purely  H-plane  (2D) 
representations.  We  then  applied  this  theory  to  guide  the  design  of  high-gain  configurations  (again, 
limited  to  2D,  H-plane  representations)  for  linear,  forward  traveling-wave,  leaky-wave  antennas.  The 
theory  built  upon  equivalent  circuit  methods  and  wave  matrix  theory,  which  provided  useful  formalisms 
upon  which  we  continue  to  build. 

During  the  second  quarter,  we  pursued  initial  extensions  of  the  previous  work  into  three  dimensions,  in 
order  to  include  phenomena  with  E-plane  dependencies.  We  succeeded  in  adding  into  the  wave-matrix 
formalism  the  refiection/transmission  properties  associated  with  the  transition  to  free  space  from  a  finite- 
width  leaky-wave  channel,  including  the  edge-tapering  essential  to  HPM  applications.  These  geometric 
aspects  do  not  arise  in  analyses  confined  to  the  H-plane  alone.  Our  3D  analyses  were  somewhat  more 
reliant  on  numerical  models  than  in  the  2D  analyses,  due  to  the  greater  complexity  of  identifying  and/or 
building  practical  analytic  approaches  capable  of  addressing  true  3D  geometries  of  interest. 

During  the  third  quarter,  we  explored  channel-to-channel  coupling  (aka,  mutual  coupling)  which  (as  we 
have  noted  earlier)  is  an  important  design  concern,  since  it  can  impact  antenna  performance  significantly 
in  terms  of  gain,  peak  power-handling,  and  impedance  matching.  Our  approach  leveraged  mostly 
numerical  methods,  along  with  some  intuitive  arguments,  as  we  explored  designs  exhibiting  different 
degrees  of  mutual  coupling  between  adjacent  channels.  As  past  and  current  antenna  literature  attest, 
mutual  coupling  analyses  are  non-trivial;  suffice  to  say,  there  is  still  much  work  to  be  done  in  this  area. 

During  the  fourth  quarter,  we  continued  to  study  and  employ  wave-matrix  based  methods,  but  with  less 
success  than  before  in  applying  this  approach  to  improve  or  optimize  the  initial  designs.  The  formalism 
itself  is  still  valid,  but  offers  reduced  practical  rewards  once  an  initial  (i.e.,  not  fully-optimized)  geometry 
(e.g.,  grill,  window,  channel  depth,  etc.)  is  derived  from  the  more  basic-level  principles.  At  that  stage,  we 
are  finding  that  further  optimization  is  currently  best  proceeding  via  numerical  means.  Additional  work 
in  the  fourth  quarter  led  us  to  identify  new  aperture  geometries  of  potentially-significant  practical  value, 
which  included  the  “BAWSEA”  and  “GAWSEA”.  These  configurations  may  significantly  extend  the 
utility  of  leaky- wave  antenna  technology  to  support  integration  on  more  challenging  platforms. 

During  the  fifth  quarter,  we  designed,  analyzed,  and  documented  representative  high-performance 
FAWSEA  and  CAWSEA  antennas  suitable  for  designation  as  “standard”  or  “recommended.”  The 
configurations  we  described  were  scalable  with  wavelength.  These  are  the  initial  entries  in  a  library  of 
antennas  that  will  continue  to  be  built  throughout  this  program. 

During  the  sixth  quarter,  we  performed  additional  investigation  of  designs  to  support  the  newer  curved 
apertures,  especially  the  “Bent  Aperture  Waveguide  Sidewall-emitting  Antenna”  (BAWSEA).  We 
presented  this  work  at  the  17th  Annual  Directed  Energy  Professional  Society  (DEPS)  Symposium  in 
Anaheim,  CA,  on  March  4th,  2015.  Our  full  slide  presentation,  entitled  “Advances  in  Low-Profile  Leaky- 
Wave  Conformable  Antennas  for  HPM  Applications,”  was  included  in  the  unclassified  proceedings  CD 
that  was  recently  distributed  by  DEPS  to  all  the  conference  attendees. 
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During  the  seventh  quarter,  we  investigated  RAWSEA  design  considerations  and  showed  that  the  angle  of 
rotation  between  the  leaky  wave  channels  and  the  aperture  can  be  understood  in  terms  of  an  equivalent 
linear  (non-rotated)  displacement,  an  interpretation  which  helps  to  guide  application  of  the  wave-matrix 
formalism.  However,  more  work  is  still  needed  to  speed-up  the  RAWSEA  design  process. 

During  the  eighth  quarter,  we  identified,  investigated,  and  applied  a  seemingly-simple  but  clarifying 
wave-mapping  methodology,  which  provided  improved  guidance  in  making  optimal  use  of  generally 
curved  platform  surfaces.  Following  this  process  helps  guide  the  designer  toward  a  solution  that  provides 
both  higher  gain  and  greater  peak  power  handling.  Via  this  approach  we  identified  and  reported  a  notable 
success  with  the  design  of  an  improved  CAWSEA  that  can  deliver  superior  gain,  yet  still  conform  to  the 
same  radius  cylinder  as  our  earlier-suggested  “standard/recommended”  design. 

For  more  information,  we  encourage  the  reader  to  refer  our  earlier  Quarterly  Reports  #1  thru  #8. 

1 .2.  Overview  of  Recent  Activities  (9th  Quarter) 

Most  of  the  technical  work  for  the  9th  quarter  was  performed  by  SARA  Staff  Scientist  Dr.  Sammuel  M. 
Jalali.  In  particular,  the  PI  (Dr.  Koslover)  asked  Dr.  Jalali  to  investigate  and  document  a  theoretical 
framework/toolset  to  help  guide  the  AAWSEA1  design  process,  in  analogy  to  the  theory  and  tools  we 
documented  in  earlier  reports  that  already  guide  the  design  of  new  FAWSEAs  and  CAWSEAs.  To  that 
end,  Dr.  Jalali  sought  to  expand  the  applicability  of  plane-wave  decomposition  (see  our  1st  Quarterly 
Report)  into  a  more  general,  but  also  inherently  more-approximate,  ray-tracing  technique.  This  approach 
employs  an  analytic  parameterization  of  the  inner-curve  (channel  back-wall)  and  outer-curve  (vicinity  of 
the  leaky-grill  wall)  ogives,  while  tracking  the  varying  angles  of  reflection  sequentially  along  the  perspec¬ 
tive  leaky  guide,  and  ultimately  adjusting  these  curves  to  yield  the  desired  output  beam.  In  principle,  this 
kind  of  approach  should  work,  given  sufficiently-gradual  curvature.  At  this  stage,  it  at  least  offers  some 
insight.  More  work  is  needed  to  advance  this  analysis,  to  integrate  it  into  the  framework  we  previously 
reported,  and  to  establish  whether  it  is  actually  useful  in  the  generation  of  practical  design  recipes  for 
efficient,  high-gain,  high-Ppk  capable,  AAWSEA  configurations. 

Further  information  about  the  aforementioned  research  is  provided  in  Section  3. 

1.3.  Upcoming  conferences 

The  DEPS  18th  Annual  Directed  Energy  Symposium,  will  be  held  in  Albuquerque,  NM,  March  7-11, 

2016.  The  Call  For  Papers  is  posted  at  http://www.deps.org/DEPSpages/DEsymp  1 6.html  with  abstracts 
due  Jan  8,  2016.  We  would  like  to  attend  this  important  conference  and  present  an  update  to  the  DEW 
community,  summarizing  the  current  status  of  the  research  and  progress  made  since  last  year.  Since  the 
abstract  deadline  is  now  rapidly  approaching,  we  will  prepare  an  abstract  for  review  by  ONR  very  shortly. 

In  addition,  we  note  that  there  is  also  an  IEEE  AP-S/USNC-URSI  conference  to  be  held  June  25-July  1, 
2016,  in  Puerto  Rico  (see  http://www.2016apsursi.org/CallForPapers.asp).  The  paper  submission 
deadline  for  that  conference  is  also  relatively  soon  (Jan.  18,  2016).  We  solicit  ONR’s  guidance  in  regard 
to  the  value  of  our  attending  and  presenting  this  work  to  the  wider-community  at  the  IEEE  AP-S/USNC- 
URSI  conference. 


1  Recall  that  the  AAWSEA,  or  “Arched  Aperture  Waveguide  Sidewall-Emitting  Antenna,”  specifically  denotes  a 
configuration  with  aperture-curvature  in  the  H-plane.  A  list  of  acronyms  we  have  used  to  date  for  variously  flat  and 
curved-surface  leaky-wave  antennas  may  be  found  in  Table  1  of  the  4th  Quarterly  Report ,  Sept.  18,  2014. 
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2.  STATUS  OF  THE  PLAN/SCHEDULE  AND  FUNDING 

Figure  1  (next  page)  maps  out  the  updated  program  plan,  for  quick  reference.  The  subject  contract  was 
awarded  on  9/18/2013  and  has  an  end  date  of  10/17/2016.  The  total  contract  value  is  $868,350,  with 
current  (per  P00005  signed  on  3/17/2015)  authorized  funding  of  $780,473.  According  to  SARA’s 
accounting  system,  as  of  Dec.  11,  2015,  expenses  and  commitments  (including  fee)  totaled  $595,502,  thus 
leaving  $184,971  available,  as  of  that  date.  If  one  simply  compares  the  calendar  and  spending  on  this 
project,  we  have  consumed  71.8%  of  the  calendar,  68.6%  of  the  total  contract  value,  and  76.3%  of  the 
authorized  funds.  We  thank  ONR  for  the  continued  support  of  this  project. 

As  the  observant  reader  may  have  noted,  some  of  the  R&D  paths  in  this  program  have  shifted  as  we  have 
gradually  identified  and  pursued  many  interesting  antenna  questions  and  novel  design  opportunities 
within  the  family  of  low-profile,  high-power-capable,  forward-traveling,  leaky-wave  antennas.  This  has 
led  to  some  schedule  slippage  of  the  planned  optimized-design  milestones  in  Task  3,  so  we  will  now  need 
to  focus  more  closely  upon  completing  those  tasks.  Overall,  there  are  no  significant  technical,  schedule, 
or  funding-related  program  problems  to  report  at  this  time. 
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Plan  of  Action  and  Milestones  (POA&M)  (updated  Dec  2015) 
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Figure  1.  Updated  Program  Plan 
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3.  RESEARCH  AND  ACTIVITIES  PERFORMED  THIS  PERIOD 

3.1 .  Generalization  of  plane-wave  wave-decomposition  to  AAWSEA  analyses 

We  consider  the  basic  analysis  of  an  AAWSEA  below.  This  configuration  is  an  important  one, 
representative  of  constraints  of  payloads  and  the  common  necessity  of  conformal  matching  of  the  antenna 
aperture  to  a  potential  DEW  platform.  The  first  part  below  is  intended  to  investigate  the  possibility  of 
mathematical  analysis  and  development  of  an  algorithm  to  compensate  for  the  phase  velocity  changes  and 
radiated  plane  wave  disturbance  caused  by  the  arched  (i.e.,  curved  in  the  H-plane)  structure.  A  simplified 
model  with  no  dielectric  window  is  considered  throughout.  In  leaky-wave  antenna  (LWA)  structures  with 
flat  side  walls  (whether  leaky  or  conducting  walls)  such  as  the  FAWSEA,  plane  waves  are  radiated  at  an 

angle  relative  to  end-fire  approximately  equal  to  the  complement  of  9a  =  cos  1  0^  with  fc  being  the 

effective  cut-off  frequency  of  the  guide  and  /  being  the  operating  frequency.  In  such  structures,  the 

guided  phase  velocity  is  constant  along  the  device  through  the  relation  of  vp  —  c,  where  Xg  is  the 

guided  wavelength  and  A0  is  the  free-space  wavelength.  For  the  purpose  of  analysis  the  guided  wave  can 
be  decomposed  in  two  interfering  TEM  plane-waves  reflecting  from  side  walls.  These  TEM  waves 
provide  a  convenient  tool  for  demonstration  and  analysis  of  the  propagating  waves.  One  of  these  waves  is 
shown  in  Figure  2  in  addition  to  the  geometric  relations  between  free-space  and  phase  velocities. 


Figure  2.  Plane-wave  representation  of  waveguide  mode. 

When  the  side  walls  are  arched,  the  phase  velocity  changes  along  the  guide  and  consequently  plane  wave 
radiation  matching  that  of  the  case  with  flat  guide  walls  no  longer  occurs.  This  result  is  consistent  with 
the  simple  plane-wave  decomposition  in  Figure  3.  In  this  figure,  both  walls  are  bent  similarly  (that  is,  the 
same  curve  is  used  for  both  side  walls).  As  can  be  seen  from  the  figure,  the  wave  direction  is  altered  and 
therefore  the  geometry  needs  to  be  modified  to  achieve  a  proper  plane-wave  output  radiation. 
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Figure  3.  Plane-wave  view  of  wave  propagation  around  an  arch  (H-plane  curve). 


3.2.  Simulating  the  guide  with  parallel  side  walls 

Let’s  first  mathematically  formulate  the  process  for  the  LWA  with  similarly  arched  side  walls  in  order  to 
demonstrate  the  existing  problem  and  verify  the  prediction  made  in  Figure  3.  Although  an  arch  with  a  true 
ogive  shape  is  desirable,  we  will  use  a  simpler  curve  equation  to  ease  our  initial  simulation.  (Ogive 
curves  are  studied  in  a  later  section.)  The  complete  simulation  is  written  in  Matlab  script  that  can  be 
simply  modified  for  other  choices  of  curves.  Let  us  assume  a  simple  quadratic  curve  for  both  side  walls 
(leaky  and  PEC),  that  is 

/(z)  =  — Ctz 2  and  g(z )  =  —  C2z2  —  h  (1) 

We  choose  a  small,  positive,  and  equal  value  for  the  quadratic  coefficients  Ct  =  C2  and  h  is  the  wave 
guide  width.  From  basic  waveguide  relations  we  have  the  cutoff  frequency  fc  =  — .  We  also  choose  an 
appropriate  operating  frequency  f  from  which  the  radiation  angle  is  computed  by 
90°  -da  =  90°  -  cos-1  (J). 

For  mathematical  analysis  we  only  include  one  of  the  plane  waves.  Obviously  the  other  TEM  wave  will 
be  treated  similarly.  Figure  4  shows  the  first  partial  reflection  of  the  TEM  wave  from  the  leaky  side  wall. 
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First  we  determine  the  equation  for  the  ray  (z)  that  passes  through  the  point  (z  —  0  ,  x  —  —h)  with 
the  slope  of 

mWi  —  tan(90°  —  9a )  =  cot£?0  (2) 

Wi(z)  =  (cot0o)z  -  h  (3) 

The  intersection  of  this  line  and  the  upper  side  wall  that  is  the  intersection  of  /  (z)  and  w 1  (z)  is 
computed  to  obtain  the  point  C^n.^n)- 

/(z)  =  ut(z)  ->  CjZ2  +  (cot0o)z  —  h  =  0  (4) 

The  roots  are 


cot0o  Vcot2  +  4 CJi 


z  =  — 


2  Ci 


2  Ci 


(5) 


cot0o  Jcot20o+4C1h 

Here  we  must  select  the  positive  root  of  ztl  = - — - 1 - — - .  Notice  that  because 

Vcot2  f?0  +  4Cx/i  >  cot£?0  the  existence  of  the  positive  root  is  guaranteed.  Having  this  intersection 
point  we  can  find  the  absolute  slope  angle  of  the  upper  side  wall. 

con  =  |tan-1(/'(zii))|  =  ltan-1(— 2C1z11)|  =  tan-1  (2  Cxzx^) 

Now  we  find  the  equation  for  the  reflected  ray  v1  (z).  The  slope  is  computed  as  follows  (see  Figure  5). 

mVi  =  —  tan(02)=  —  tan(oj11  +  (90°  —  (pi)) 

<Pi  =  90°  -  (90°  -  0O )  -  0JL1  =  eo  -  0)lt 


mVi  —  —  tan(90°  —  90  +  2£on)  —  —  cot(9a  —  2a)n)  (6) 

The  ray  (z)  also  passes  through  the  point  (z±±,xlt)  and  therefore  the  line  equation  can  be  written 
using  (6). 

z)  =  mVi(z  -  zn)  +  Xu  =  —  cot(0o  -  2£un)z  +  x±1  +  cot (9a  -  2co11)z11  (7) 
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Now  equipped  with  the  equation  of  the  ray  v1  (z)  and  the  lower  side  wall  equation  g(z)  we  find  the 
intersection  point  (z12/  x12)  shown  in  Figure  5. 

g(z)  =  vt(z)  ->  — C2z 2  —  h  =  —  cot(0o  —  2a)lt ) z  +  xtl  +  cot (90  —  2a)tl ) z1]L  (8) 


Equation  (8)  leads  to  another  quadratic  equation  that  is  being  solved  to  find  the  intersection  point 

(Zl2^12)- 

C2z 2  -  cot(0o  -  2a)11)z  +  [h  +  xtl  +  cot(dQ  -  2oj11)z11]  =  0  (9) 

Equation  (9)  has  two  roots  with  at  least  one  positive  (with  the  similar  justification  as  for  equation  (4)).  Let 
us  denote  this  root  by  z12  from  which  the  x12  —  g(z12)  —  —C2z\2  —  h  is  computed.  This  point  can  be 
used  to  find  the  absolute  slope  angle  of  the  lower  side  wall,  per  Figure  5. 

co12  =  |tan_1(^'(z12))  I  =  |tan-1(— 2 C2z12)  |  =  tan-1(2C2z12) 

Also  notice  that  the  angle  d2  is  known  from  the  slope  of  (z)  in  previous  step. 

92  =  90°  -  90  +  2 (On 

ip1  —  90°  —  ( 02  —  CJ12)  =  90°  —  02  +  o)12  —  9a  —  2c otl  +  oi  12 

The  critical  angle  that  also  specifies  the  slope  of  the  reflected  ray  w2  (z)  in  Figure  4  is  denoted  by  6 '  and 
is  found  by  the  following. 

r  =  90°  -Tp1-  co12  =  90°  -  Q0  +  2 <yn  -  2£U12  (10) 


Now  the  equation  for  the  third  reflected  ray  can  be  computed  using  the  slope  angle  9  *  and  the  intersection 
point  (z12,  x  12). 


w2(z)  -  ttv2(z  -  Z12)  +  *12 

mw,2  =  tan(0*)  =  tan(90°  —  90  4-  2cjx1  —  2<y12)  =  cot(0o  —  2a )n  4-  2cj12) 

w2(z)  =  cot(<90  -  2a)n  +  2 oj12)z  +  [x12  -  cot(da  -  2a)n  +  2 a)12)z12]  (11) 

Finally  this  process  can  be  carried  on  for  as  many  reflections  as  there  are  in  the  guide  length  not  only  for 
the  plane  wave  shown  in  the  above  but  also  for  the  complementary  plane  wave  that  starts  from  the  top 
side  wall. 


We  have  prepared  a  Matlab  script  (included  here,  for  convenient  reference)  to  perform  these  operations 
and  plot  the  implied  radiated  rays,  to  visualize  the  problem  that  arises  from  variation  of  the  phase  velocity 
in  an  arched  leaky  wave  guide.  The  parameters  used  are:  h  =  0,1651  m,  L  —  1.0  m,  and  /  =  1  GHz. 


%  Script  name:  "Vp_testl.m"  (C:\Matlab\NAV41) 

%  This  script  is  used  to  demonstrate  the  effect  of  arched  guide  on  the 
%  phase  velocity  and  on  the  radiated  plain  wave. 


%  Sam  Jalali,  SARA,  Inc.,  Cypress,  California 
%  Version:  1.0,  10/30/2015 

%  ************************************************************************* 
close  all;  clear  all; 
format  long;  format  compact; 


c_const  = 

L  =  1.0; 

h  =  1 6  5 . 1  e  -  3 ;  % 

Cl  =  0.10;  % 

C2  =  Cl;  % 

fc  =  c_const/ (2*h) ; 
f  =  1 . 0  e  9  ; 

Tets  =  acos(fc/f); 

%  Draw  the  upper  and 
z=l in space ( 0 , L, 500 ) ; 
for  n  =  l:length(z) 
XI  (n)  =  - C 1  * z  (n) 
X2  (n)  =  -C2*z (n) 


99792458e8; 

%  guide  length  in  meters 
guide  width  in  meters 
upper  side  wall  equation  f(z)  = 
lower  side  wall  equation  g(z)  = 
%  cutoff  fequency  in  GHz 
%  operating  frequency 
%  TEM  modes  angle 
lower  side  walls  of  the  guide 


-Cl*z/ 

-C2*z/ 


2 

2-h 


2; 

‘2-h; 


upper 

lower 


sidewall 

sidewall 


(grid)  equation 
equation 


end 

plot ( z , XI , ' :k', ' LineWidth ',3.5);  hold  on; 
plot ( z , X2 ,  ' -k ' ,  ' LineWidth '  ,  3 . 0 )  ;  grid  on; 

xlabel ( ' z '  ,  ' Font Size '  ,  18 ,  1 FontWeight ' ,  ' bold ' ,  'Color ' ,  ' b ' ,  ' Font Angle ' ,  .  .  . 

' italic FontName Times  New  Roman'); 
ylabel ('x', 'FontSize',18, ' FontWeight ', 'bold', 'Color', 'b', ' FontAngle ' , . . . 

' italic ',' FontName ',' Times  New  Roman'); 
title ( ' TEM  mode  decomposition ' ,  ' Font Size '  ,  14  ,  ' FontWeight ' ,  ' bold ' ,  .  .  . 

' Color '  ,  ' b '  ,  ' FontName ' ,  ' Times  New  Roman ' ) ; 


N  =  7;  %  number  of  reflections  in  the  guide 

%  Start  of  the  loop  (N  iterations) ,  for  lower  TEM  mode  ray 
zo  =  0;  Ko  =  -h;  TetO  =  Tets; 
for  n=l : N 

%  Compute  the  intersection  of  wl (z)  with  upper  side  wall 
ml  =  cot (TetO);  %  slope  of  wl (z) 

ZR  =  roots ( [Cl , ml , Ko-ml*zo] ) ; 
zl  =  Z  R ( 2 )  ; 
if  zl  <=  0 

error (' The  z  vale  is  negative'); 

end 

xl  =  -Cl*zlA2; 

%  Draw  the  ray  section  wl (z) 

XL  =  1.0;  %  top  X  limit  for  plane  wave 

ZL  =  (XL-Ko+ml*zo) /ml ;  %  corresponding  Z  limit 

if  ZL  >=  L 
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ZL  =  L; 

end 

z  =  linspace (zo, ZL, 300) ; 
for  n  =  l:length(z) 

wl (n)  =  ml*z (n) +Ko-ml*zo; 

end 

plot  ( z , wl , '  -  r '  , ' LineWidth ',3.0) ; 

%  Compute  the  next  ray  section 

Tetl  =  abs  (atan (-2*Cl*zl) ) ;  %  wl (z)  slope  angle  absolute  value 

%  Phil  =  TetO-Tetl;  Phi2  =  pi/2-Phil+Tetl ; 

Phi2  =  pi/2 -TetO+2 *Tetl ;  %  angle  of  w2(z)  in  absolute  value 

%  Compute  the  intersection  of  w2(z)  with  lower  side  wall 
m2  =  -tan(Phi2);  %  slope  of  w2(z) 

K1  =  h+xl-m2*zl; 

ZR  =  roots ( [C2,m2, Kl] ) ; 
z2  =  ZR ( 2 ) ; 
if  (z2  <=  0) 

error ('The  z  vale  is  negative'); 

end 

x2  =  -C2*z2^2-h; 

%  Draw  the  ray  section  w2(z) 
z  =  linspace ( zl , z2 , 300 ) ; 
for  n  =  l:length(z) 

w2 (n)  =  m2*z (n) +xl-m2*zl ;  %  (z2,x2)  could've  been  used  too 

end 

plot  ( z , w2 ,  ' -r '  ,  ' LineWidth ',3.0)  ; 

%  Prepare  the  parameters  for  next  iteration 

Tet2  =  abs  (atan (-2*C2*z2 )) ;  %  w2(z)  slope  angle  absolute  value 

%  Psi  =  pi/2-Phi2+Tet2 ;  Tets  =  pi/2-Psi-Tet2 ; 

Tetc  =  pi/2-TetO+2*Tetl-2*Tet2 ;  %  the  crtical  angle 
TetO  =  pi/2-Tetc;  %  New  TEM  angle 

zo  =  z2;  Ko  =  x2; 

end 

%  Start  of  the  loop  (N  iterations),  for  upper  TEM  mode  ray 
zo  =  0;  Ko  =  0;  TetO  =  Tets; 
for  n=l:N 

%  Compute  the  intersection  of  wl (z)  with  upper  side  wall 
ml  =  -cot (TetO);  %  slope  of  wl (z) 

ZR  =  roots ( [C2, ml, Ko-ml*zo+h] ) ; 
zl  =  ZR ( 2 ) ; 
if  zl  <=  0 

error (' The  z  vale  is  negative'); 

end 

xl  =  -C2*zl"2-h; 

%  Draw  the  ray  section  wl (z) 
z  =  linspace (zo, zl, 300) ; 
for  n  =  l:length(z) 

wl (n)  =  ml*z (n) +Ko-ml*zo; 

end 

plot ( z , wl , ' -m ' , ' LineWidth ',3.0) ; 

%  Compute  the  next  ray  section 
Tetl  =  abs (atan (-2*C2*zl )) ; 

%  Phil  =  pi/2-TetO;  Phi2  =  pi/2-Phil+Tetl ; 

Tetc  =  pi/2 -TetO-2 *Tetl ;  %  angle  of  w2(z)  in  absolute  value,  criitical 

%  Compute  the  intersection  of  w2(z)  with  lower  side  wall 
m2  =  tan (Tetc) ; 

Kl  =  xl-m2*zl ; 

ZR  =  roots ( [Cl, m2, Kl] ) ; 
z2  =  ZR  (2  )  ; 
if  (z2  <=  0) 

error (' The  z  vale  is  negative'); 

end 

x2  =  -Cl*z2^2; 

%  Draw  the  ray  section  w2(z) 

XL  —  1.0;  %  top  X  limit  for  plane  wave 

ZL  =  (XL-Kl)/m2;  %  corresponding  Z  limit 

if  ZL  >=  L 
ZL  =  L; 

end 

z  =  linspace ( zl , ZL, 300 ) ; 
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for  n  =  l:length(z) 

w2 (n)  =  m2*z (n) +xl-m2*zl ;  %  (z2,x2)  could've  been  used  too 

end 

plot  ( z , w2 ,  ' -m ' ,  ' LineWidth ',3.0) ; 

%  Prepare  the  parameters  for  next  iteration 

Tet2  =  abs (atan (-2*Cl*z2 ) ) ;  %  wl (z)  slope  angle  absolute  value 

TetO  =  pi/2-Tetc-2*Tet2;  %  New  TEM  angle 

zo  =  z2;  Ko  =  x2; 

end 


Figure  6  shows  the  result  of  executing  the  above  Matlab  script.  As  can  be  seen,  some  compensating 
technique  is  required  to  return  the  rays  to  parallel.  And  that  means  an  adjustment  to  the  phase  velocity,  as 
expected,  is  necessary. 


Figure  6.  Departure  from  parallel-ray  outputs  due  to  the  curved  walls. 


3.3.  Finding  the  side  wall  equation  for  proper  phase  velocity  compensation 

The  above  procedure  was  based  on  two  simplifying  assumptions.  First  the  leaky  wall  shape  was 
considered  to  be  a  simple  quadratic,  which  for  practical  reasons  might  be  changed  to  ogive  or  Bezier 
shapes.  Second,  the  lower  PEC  wall  was  assumed  to  be  known  and  parallel  to  the  leaky  wall  (same 
equation  with  an  offset).  In  fact,  the  lower  PEC  wall  should  actually  be  computed  for  the  proper 
compensation. 

Here  we  retain  the  same  simple  curvature  for  both  the  leaky  wall  and  the  PEC  wall  to  develop  the 
methodology  in  a  simple  way.  Generalization  of  the  methodology  will  be  discussed  later.  Therefore,  once 
again  we  assume  the  leaky  wall  to  have  the  equation  / (z)  =  —  Ctz2  and  the  lower  side  wall  have  the 
same  simple  quadratic  equation  g(z)  =  —  C2z2  —  h.  However  this  time  C2  =£  Ct  and  in  fact  we  need  to 
compute  the  appropriate  value  for  C2  (if  possible). 
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Refer  again  to  Figure  5.  For  proper  operation  we  require  that  the  reflected  ray  w2  ( z )  be  in  parallel  with 
the  first  ray  in  Figure  4  that  is  w^z)  and  similarly  towards  the  end  of  the  guide.  Such  requirement 
mandates  that  the  slope  angle  of  w2(z)  that  is  8*  be  equal  to  the  slope  angle  of  w1(z)  that  is  90°  —  80. 
In  fact,  this  was  the  reason  we  called  9  *  the  critical  angle. 

Now  we  will  proceed  to  show  how  the  previous  requirement  leads  to  a  method  of  finding  the  proper  value 
of  C2  that  determines  the  shape  of  lower  compensating  wall.  We  have  the  critical  angle  value  from 
9  *  =  90°  —  9Q.  From  equation  (10)  in  previous  section  we  proceed  to  find  the  slope  of  the  tangent  line 
to  the  lower  PEC  wall  a)12  at  the  point  of  intersection  (z12,  x12)  that  is  yet  to  be  determined  in  this  case. 

8 *  =  90°  -eo  =  90°  -90  +  2a)lt  -  2a)12  (10) 

So  the  value  for  this  angle  can  be  computed  to  be  equal  to  the  known  angle  <un  in  the  following  equation 

a)12  =  a)lt  (12) 

We  also  know  from  the  slope  of  the  tangent  line  to  the  lower  side  wall  with  g{z)  =  —  C2z2  —  h  that 
co12  =  ||tan_1(g'(z12))  |  =  |tan_1(-2C2z12)  |  =  con  (13) 

From  the  equation  (13)  we  can  compute  the  intersection  point  z12  in  terms  of  the  unknown  coefficient  C2 
by  taking  “tan”  of  the  both  sides  as  follows. 

tan(aji  t ) 

2C2z12  =  tan(cjn)  ->  z12  =  — — -  (14) 

But  this  value  of  z12  must  satisfy  the  equation  (9)  as  follows. 

C2z\2  —  cot (80  —  2 a)lt)  z12  +  [h  +  xlt  +  cot (90  —  2a>tl)  zn]  =  0  (15) 

Inserting  equation  (14)  into  (15)  we  can  compute  C2  as  follows. 

c  _  tan(ojn)  [2  cot(fl0  -  2ajn)  -  tan(qjn)] 

4[^1X  +  h  +  cot(6a  -  2to11)z11] 

The  above  procedure  demonstrate  a  general  methodology  used  to  determine  the  geometrical  aspects  of  the 
lower  wall  of  the  LWA  to  compensate  for  phase  velocity  variation  in  AAWSEA.  It  also  indicates  the 
complications  involved  in  the  design  of  an  arched  aperture  vs.  a  flat  one. 

In  the  general  case,  one  can  use  this  technique  to  find  several  computed  points  for  the  lower  PEC  side 
wall.  Although  the  number  of  computed  points  is  limited,  an  interpolation  of  these  points  gives  an 
equation  in  terms  of  polynomial  that  can  be  a  good  approximation  for  the  target  compensating  curve. 

A  modified  Matlab  script  (Vp_Test2.m)  is  used  to  demonstrate  the  result  of  this  methodology.  Figure  7 
shows  the  correction  made  only  by  computing  one  compensation  point  on  the  lower  PEC  wall.  Obviously, 
the  compensation  must  be  performed  for  many  points  (or  more  generally,  a  continuum)  along  the  guide. 
An  appropriate  analytic  function  or  a  high-order  polynomial  can  likely  be  employed  for  satisfactory 
results.  We  will  discuss  the  latter  in  the  next  section. 
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3.4.  Ogive  curvature  and  practical  restrictions 

Use  of  ogive  curves  is  common  in  minimizing  the  drag  caused  by  air  resistance.  We  will  consider  the 
tangent  ogive  curve  given  by  the  following  mathematical  formula  written  for  our  previously  used 
coordinate  system. 


/O) 


In  this  formula,  d  is  the  curvature  diameter  and  C  —  -  is  called  the  caliber.  Such  an  ogive  curve  has  an 

offset  of  / (0)  =  In  order  to  match  the  curve  with  our  upper  side  wall  with  no  offset  of  leaky  guide 
(/ (0)  =  0)  we  modify  and  simplify  the  formula  to  the  following: 


/(z)  = 


a  (17) 


The  main  parameter  in  the  above  ogive  curve  is  a  —  d  yC2  +  -J.  The  methodology  for  compensating 

phase  velocity  via  the  geometric  optics  technique  (used  in  previous  sections)  needs  to  be  performed 
iteratively  for  multiple  reflections.  Figure  8  is  an  extension  of  Figure  4  and  Figure  5.  It  is  used  here  to 
clarify  the  detailed  procedure. 
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Figure  8.  Extending  the  discrete  ray  tracing  process. 

Similar  geometrical  derivations  are  applied,  except  that  now  the  leaky  guide  wall  has  an  ogive  curve  and 
we  also  assume  that  the  lower  PEC  wall  can  have  an  ogive  curve  with  different  parameter,  namely: 

g(z)  —  V b 2  —  z2  —  b  —  h  (18) 

Notice  that  for  larger  diameter  (lower  curvature)  of  the  PEC  wall  we  need  to  have  b  >  a,  moreover,  with 
very  high  parameter  value  ( b  »  1)  the  curve  becomes  almost  horizontal. 

The  similar  derivation  for  shape  function  is  summarized  in  the  following. 

a.  Determine  the  intersection  of  the  ray  specified  by  w^z)  and  the  leaky  wall  f(z ). 

w1(z)  =  tan  (90°  —  dQ)z  —  h  =  cot  (60)z  —  h 
f(z)  =  V  a2  —  z2  —  a 

Solve  for  the  nonlinear  equation  V&2  —  z2  —  a  —  cot(0o)  z  +  h  =  0,  which  can  be  done  by 
“fzero”  function  in  Matlab  to  obtain  zt  1 .  And  the  corresponding  coordinate 

xlt  =  cot (80)  zlt  —  h 

b.  Compute  the  absolute  of  slope  angle  of  f  (z)  at  the  intersection  point. 

/'(z)  =  -  =  ,  GJn  =  tan-1(— /'(zii)) 

V  a  ~  zii 
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c.  Find  the  equation  of  the  reflected  ray  V\(z)  as  follows. 


<Pi  =  90°  -  (90°  -  eo)  -  mn  =0O-  (olt 
mVi  =  —  tan(mn  +  90°  —  (p{)  =  —  tan(90°  —  0O  +  2(o11')  =  —  cot(0o  —  2a in) 
vx(  z)  —  —  cot(0o  —  2a>n)z  +  x11+  cot(0o  —  2a)11)z11  (19) 

d.  The  intersection  point  (z12,  *12)  of  the  latter  ray  v1  (z)  and  the  lower  side  wall  specified  by 
g(z)  cannot  be  done  because  of  the  unknown  parameter  b  in  the  lower  side  wall  equation.  If 
we  had  this  intersection  point  the  new  absolute  slope  angle  could  be  computed  by 

cj12  =  |tan_1(y(z12))|  -  tan_1(— 5'(z12))  -  tan-1  (q==  2  j  (20) 

But  according  to  the  same  process  that  culminated  in  equation  (13)  we  should  require  that  the  slope 
angle  be  equal  such  that  the  leaky  wave  be  in  parallel.  That  is  a)12  =  Cd1:L. 

This  requirement,  in  conjunction  to  equation  (20)  for  known  Ct)ll9  can  be  used  to  compute  the 
abscissa  of  intersection  point  z12 . 

tan(“ll) = w^w2  <21) 

Solving  the  equation  (21)  for  z12  with  known  value  of  a)lt  we  will  have 

z12  =bsm(a)tl )  (22) 

Obviously  z12  =  b  sin(6)lx)  >  0  is  guaranteed.  The  corresponding  function  value  of  this 
intersection  point  can  also  be  found  (after  some  algebra)  in  terms  of  unknown  b  as 

*12  =  9 (^12)  =  g(b  sin(con))  =  ij(cos(con)  -  1)  (23) 

Here  we  also  see  that  x12  =  j&tcostcd^)  —  1)  <  0  because  costo^)  —  1  <  0  as  we  require  for 
proper  value  of  x  ordinate. 


tUi2  =  CL>n  =  tan  1 


z12 


2  —  72 

z12' 


e.  Finally,  the  intersection  point  (z12, x12)  that  are  computed  in  terms  of  b  value  must  satisfy  the 
equation  of  v1  (z)  in  (19). 


*12  =  ^l(z12)  =  —  cot(£?0  -  2cj11)z12  +xL1  +  cot(eo  -  2(0 11)z11 
^(costcun)  —  1)  =  —  cot(f?D  —  2 (0^)  ( b sintaq!))  +  xxl  +  cot(0o  —  2cjn)  zlx 


Solving  for  unknown  b  value  we  have 

xri  +  cot (9a  -  2(on)  +  zn 
b  — - - - - - 

COstcd-Li)  +  sin^!!)  cot(0o  —  2con)  —  1 


(24) 


The  computed  value  of  b  for  this  intersection  point  determines  the  main  parameter  of  the  ogive 
shape  we  approximate  for  the  lower  wall,  or  at  least  provide  us  with  a  point  (z12,  x12)  on  that  wall 
that  can  be  collected  with  further  intersection  point  values. 


f.  The  process  outlined  in  the  above  steps  must  be  repeated  to  obtain  many  points  of  intersection 
(zi2,*£2)  for  i  —  1, 2, 3, ...  as  shown  in  Figure  8  for  two  iterations.  The  collection  of  these 
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points  can  be  used  to  determine  the  shape  of  the  lower  PEC  side  wall  of  the  leaky  guide.  We  will 
use  a  Matlab  script  to  demonstrate  this  process. 


g.  Based  on  the  derivations  of  the  above  steps,  we  can  recognize  a  theoretical  limit  for  the  bend 
curvature  depth  beyond  which  the  wave  stops  in  the  guide  and  behaves  similar  to  the  cutoff 
frequency.  The  slope  of  the  partially  reflected  wave  i^(z)  from  the  leaky  side  must  remain  less 
than  90°  in  absolute  value  for  the  TEM  wave  to  move  forward  in  the  guide  that  is 
90°  —  9Q  +  2o)11  <  90°  (see  equation  (6)  or  the  one  before  (19)). 


— 90  "h  2(x)^  -jl  <  0  — >  co <-  —— 

2 


Or  generally  for  each  reflection  point 


-0o  +  2(x>il  <  0  ->  a)it  < 


9, 


(25) 


But  we  know  that  0o  —  cos 


"(7) 


and  for  any  slope  angle  in  the  iterative  computation  we  have 

60^  =  tan-1  (  —  /'(Zii))  with  i  being  the  index  of  the  iteration  step.  Then  using  these  relations  we 
have 


cos  1  (y)  >  2 tan  ‘(-f'(2;,))  =  2tan  1  f  = 
f 


Taking  the  cosine  of  the  above  equation,  we  have 


fc 

7 


>  cos  2  tan 


-i 


z;i 


W 


-2  -  4 


!) 


(26) 


Now  using  the  identity  of  cos (2 a)  =  1+tan2  ^  we  Proceed  with  the  equation  (26)  as  follows. 


cos  2  tan 


-l 


Zil 


\\ 


1  - 


zil 


a 2  —  zfi  a2  -  2  zf 


il 


\Ja2-4// 


1  + 


(27) 


az  -  z, 


il 


The  equation  (26)  and  (27)  gives  us  a  limitation  on  the  operating  frequency  in  relation  to  the  size  or 
depth  of  curvature. 


fc  a  -  2z/ 

t  ^ 


a* 


f<fc 


a 


fi2  -  zL 


(28) 


Now  if  we  use  the  maximum  desirable  value  for  z  that  is  the  length  of  leaky  guide  L  and  replace  the 


ogive  parameter  (and  the  formula  for  the  caliber)  a  =  d  ^C2  +  ^ 
frequency  limitation  in  terms  of  the  ogive  length  and  diameter. 


4d 


■  we  can  translate  the 
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f  <fc 


(4  ]}  +  d2) 


2\2 


f<fc 


( 4L2  +  d2\2 

(29) 


/4 L2  +  d2)2  -  1 6L2d2/ 

The  following  Matlab  script  performs  this  compensation  technique  for  the  ogive-shape  wall. 


%  Script  name:  MVp_test3.m" 

%  This  script  is  used  to  demonstrate  the  methodology  of  finding  the  points 
%  on  the  lower  side  wall  when  the  upper  leaky  wall  id  an  ogive  with  the 
%  parameters  in  the  script  and  the  lower  wall  is  approximated  similarly. 


%  Sam  Jalali,  SARA,  Inc.,  Cypress,  California 
%  Version:  1.0,  11/05/2015 

%  ************************************************************************* 
close  all;  clear  all; 
format  long;  format  compact; 


L  =  1.0;  %  cone  length 

d  =  3.0;  %  cone  base  diameter 

C  =  L/d;  %  caliber 

a  =  d* (CA2+l/4)  ; 
c_const  =  2 . 997 92458e8 ; 

h  =  165.1e-3;  %  guide  width  in  meters 

fc  =  c_const/ (2*h) ;  %  cutoff  fequency  in  GHz 

f  =  1 . 3e9 ;  %  operating  frequency 

Teto  =  acos (fc/ f ) ;  %  TEM  modes  angle 

ml  =  cot (Teto) ; 


%  Plot  the  complete  upper  side  wall 
z  =  linspace  (0,L,200); 
for  n  m  1:200 

x(n)  =  sqrt (aA2-z (n) A2 ) -a; 

end 

plot ( z , x, ' : k ' , 1 LineWidth ' , 3 . 0 ) ;  hold  on; 
mZ  =  0;  X  =  -h;  Z  =  0; 

N  =  10;  %  number  of  reflections/iterations 

zl  (1)  =  0;  xl  (1)  =  0; 
z2  (1)  =  0;  x2  (1)  =  -h; 
wl (1)  =  0; 

%  main  loop 
for  n  =  2:N 

%  W2 (n)  =  ml*z+x2 (n) -ml*z2 (n) 

%  Compute  the  first  intersection  point,  zll,  xll 
funl  =  @  (z)  sqrt  (a^-z^C  )  -ml*z-a-X+mZ ; 

zl (n)  =  f zero (funl, [Z, L] ) ; 
xl (n)  =  ml*zl (n) +X-mZ ; 

%  xl (n)  =  sqrt (aA2-zl (n) ^2 ) -a;  %  alternative  equation 

fpzl  =  -zl (n) / sqrt (aA2-zl(n)A2);  %  drivative  value,  f '  (zl  (n) ) 

wll  =  atan(-fpzl); 

wl (n)  =  wll;  %  slope  angle  is  stored  to  testing 

%  Stop  criterion  check 
if  wll  >=  Teto/2 
break; 

end 

%  Compute  the  second  intersection  point,  zl2,  xl2 

Cl  =  cot (Teto-2*wll ) ;  C2  =  cos (wll);  C3  =sin(wll); 

b  (n)  =  (xl (n) +Cl*zl (n) +h) / (C2+C3*C1-1) ; 

z2 (n)  =  b (n) *sin  (wll ) ; 

x2  (n)  =  b  (n)  *  (cos  (wll) -1) -h; 

mZ  =  ml*z2  (n)  ;  X  =  x2  (n)  ;  Z  =  z2  (n)  ; 

end 

%  plot ( zl , xl ,': ms LineWidth ', 3 ,' MarkerSize ', 4 ) ;  hold  on; 
plot ( z2 , x2 ,' -ms LineWidth ', 3 ,' MarkerSize ’, 4 ) ;  grid  on; 

for  n  =  1 : 4 

line ( [ z2 (n) , zl (n+1 ) ] , [x2 (n) , xl (n+1 ) ] ) ; 
line ( [ zl (n+1 ) , z2 (n+1 ) ] , [xl (n+1 ) , x2 (n+1 ) ] ) ; 

end 
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The  result  of  executing  the  above  script,  with  some  additions  and  observing  the  limitation  discussed  in 
part  “g”  above,  is  shown  in  Figure  9.  Here  the  curvature  is  somewhat  exaggerated  for  demonstration 
purposes.  Practical  curves  are  significantly  smoother.  The  discussion  in  part  “g”  explains  why  the  guided 
wave  is  not  carried  on  past  a  certain  z  point.  Appropriate  parameters  can  be  used  for  complete  guide 
length  simulation. 


3.5.  Design  approximation  for  the  lower  (PEC)  wall  of  the  LWA 

In  this  section  we  attempt  to  develop  and  test  an  algorithm  for  computing  grid  wire  diameters  for  an 
AAWSEA  similar  to  what  has  been  successfully  done  for  the  FAWSEA.  Since  the  geometry  of  the 
AAWSEA  makes  the  required  steps  for  this  part  of  project  more  complicated,  we  do  not  include  a 
dielectric  window  initially. 

The  leaky  wall  of  the  assumed  guide  will  be  curved  according  to  ogive  shapes  with  zero  offset  as 
described  in  Section  3.3.  Generally,  for  waveguide  curvatures  where  the  radius  is  greater  than  20  times 
the  wavelength,  variations  of  important  parameters  such  as  wave-number  are  slow,  and  some  appropriate 
approximations  are  justified.  We  will  take  advantage  of  that  here,  to  simplify  the  algorithm,  while  leaving 
detailed  verification  to  numerical  models. 

A  mathematical  analysis  was  proposed  in  section  3.4  in  order  to  compensate  for  the  AAWSEA  curvature, 
such  that  an  almost  plane  radiated  wave  can  be  achieved.  As  a  result,  several  points  on  the  lower  PEC 
wall  were  considered.  A  study  of  these  points  did  not  reveal  a  simple  analytical  functional  representation. 
Various  interpolating  polynomial  functions  are  available  through  the  standard  numerical  methods; 
unfortunately  employing  such  a  polynomial  resulted  in  significant  error  in  the  slope  of  the  tangent  lines  at 
the  interpolating  points. 
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Regarding  the  quest  for  a  more  appropriate  approximation  to  the  curve  of  the  compensating  PEC  wall,  it 
should  be  noted  that  these  ogive  curves  can  fit  only  one  point  perfectly  (aside  from  the  common  starting 
point).  The  most  significant  errors  occur  towards  the  end  of  the  LWA  (see  Figure  10). 


z-axis 

Figure  10.  Ogive  curves  with  different  parameters,  with  a  common  starting  point. 


Three  different  alternatives  were  tested:  A  function  created  by  interpolating  a  polynomial  that  passes 
through  all  the  computed  points  seemed  to  be  a  natural  option,  however  this  function  will  have  significant 
error  in  tangent  line  slopes  at  the  interpolating  points.  The  second  option  was  an  ogive  curve  with 
matching  parameters  for  the  end  point  of  the  LWA  and  the  third  was  a  simple  quadratic  function  also 
compensating  for  the  end  point.  A  Matlab  script  was  written  to  test  the  performance  of  these  options.  The 
script  finds  the  reflection  angles  by  solving  for  the  intersection  points  as  outlined  in  the  previous  section 
and  shown  in  Figure  8.  Studying  the  ogive  approximation  versus  the  quadratic  one  showed  that  the 
quadratic  function  performs  better,  with  a  maximum  reflection  angle  error  of  1.24  deg  and  minimum 
angle  error  of  0.23  deg,  when  the  ogive  curve  fit  had  a  maximum  reflection  angle  error  of  1.21  deg  and 
minimum  of  0.49  deg. 

Figure  1 1  shows  how  closely  the  quadratic-fit  curve  for  the  lower  PEC  wall  follows  the  computed  points 
(the  points  are  also  shown).  Here  the  leaky  wall  ogive  function  was  selected  for  a  guide  of  about  1  meter 
long,  cone  base  diameter  of  14  meters  with  the  following  formula.  (Figure  1 1  indicates  a  short  flat  input 
section  guide  as  well.) 

/(z)  =  V«2  -z2  -  CL,  a  —  3.5714286  (30) 

The  compensation  achieved  by  the  quadratic  approximation  g{z )  =  CQz 2  —  h  can  be  seen  in  Figure  12 
where  the  plane  wave  radiation  is  predicted.  As  a  conclusion,  we  simplify  the  complete  analysis  of 
AAWSEA  used  for  determining  the  grid  wire  diameters  by  assuming  perfect  compensation  for  reflection 
angle  and  similarly  for  phase  velocity.  The  latter  is  an  important  assumption,  which  not  only  eliminates 
complicated  angle  computations  as  outlined  in  the  previous  section  but  also  will  be  applied  in  some  other 
parameter  calculations  as  will  be  explained  in  the  following  sections.  Symbolically  we  refer  to  this 
assumption  by  di  =  0O  for  i  =  1,  2, ... 
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Figure  11.  Use  of  a  quadratic  fit  curve  for  the  back  wall. 


Figure  12.  Predicted  rays,  radiated  plane  wave. 
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3.6.  Mapping  from  flat  axis  to  the  arched  line 

A  FAWSEA  possesses  a  flat  leaky  wall  along  the  z-axis.  After  introducing  an  arch,  according  to  an 
ogive,  we  have  the  configuration  known  as  an  AAWSEA  (Figure  13). 


Figure  13.  Predicted  rays,  arched  configuration 


Let  us  calculate  the  arc  length  denoted  by  ^  =  /  (z)  on  the  arched  wall  specified  by  the  equation  (30).  By 
doing  this  calculation  each  point  on  the  z-axis  can  be  mapped  to  a  corresponding  point  on  the  arched  wall. 


fTd<  =  /„Vdz2  +  d*2T  J1  +  (s)  dz  =  /0'/TTF 


f'(z)  =  - 


Va2  — 


\dzj 

i  +  [TOO]2  = 


'(z)]2  dz 


a 


a2  —  z2 


£  —  l  - d —  a  sin  1  (— i  =  a  sin  1  (— )  (31) 

Jo  ^a2  -  Va/  o  Va^ 

With  the  conversion  from  the  z-axis  to  the  arc  length,  the  reverse  conversion  is  also  known. 


z  —  a  sin 


m(-)  (32) 


3.6.1 .  Wave  leakage  rate  on  the  arched  wire  grid 

An  ideal  attenuation  rate  that  represents  the  wave  leakage  from  the  guide  was  formulated  for  FAWSEA 
previously  (e.g.,  see  the  1st  Quarterly  Report).  This  attenuation  rate,  which  varies  with  distance  along  the 
guide,  has  been  proved  to  be  successful  and  can  be  generalized  for  AAWSEA.  The  ideal  attenuation  and 
its  conversion  to  the  arc -variable  are  given  in  the  following  equation. 

1  _  _  1 

7  ~  r  ^ideal  CO  —  TFT  (^3) 

L-asin(|) 
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3.7.  Ideal  power  reflection  coefficient 

Another  important  formalism  successfully  tested  in  analysis  of  FAWSEA  antenna  is  related  to  derivations 
of  the  ideal  power  reflection  along  the  guide  based  on  the  first  order  differential  equation  formulation. 

Rp, Ideal  00  =  exp (~Saideal(z))  (34) 


In  the  above  equation  S  represents  the  distance  between  consecutive  intersections  of  plane  wave  rays  on 
the  grid,  and  is  constant  for  FAWSEA. 


S 


o 


4ft2  4ft2  sin  eo 


(35) 


In  equation  (35)  the  90  is  the  complement  of  the  desired  radiation  angle  and  h  is  the  guide  width,  per 
Figure  14. 


This  parameter  h  is  no  longer  constant,  in  contrast  to  the  flat  aperture  case  (FAWSEA).  It  will  vary  from 
point  to  point  along  the  guide  as  seen  in  Figure  13.  First,  notice  that  we  require  a  minimum  of  flat  input 
section  (no  leakage)  for  the  guide  that  is  determined  by  £0  =  h  tan  90  for  proper  injection  of  the  wave 
into  the  leaky  guide.  Now  in  order  to  achieve  a  good  estimate  of  this  length  denoted  by  S(  we  can 
compute  the  successive  intersecting  points  of  the  rays  inside  the  guide  as  shown  by  points  (zilt  xit)  and 
(zi+l  lrxi+11),  and  the  approximation  for  5)  is  done  by  computing  the  straight  distance  between  these 
points,  as  in: 

f  2  2 

Si  =  ^(zi+i,i  ~  z;i)  +  ~  Xu)  (36) 

The  concise  methodology  used  to  compute  these  points  is  summarized  in  the  following. 

•  The  equation  for  the  leaky  wall  f  (z)  =  Va2  —  z2  —  a  and  the  PEC  (lower)  wall 
g(z )  =  Caz2  —  h  are  known. 
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•  We  will  place  the  grid  wires  at  equal  distance  on  the  arched  wall  denoted  by  Q  for  i  =  1,2, ... 

•  For  each  location  on  the  arched  wall  ^  we  compute  the  starting  point  (zi+l  l,xi+l  l). 

•  Using  the  equation  for  the  plane  wave  ray  wi+1(z)  =  —  zi+1  x)  +  xi+11  when 

mw  =  cot 9i  =  cot6a  and  g(z )  =  Caz2  —  h  is  known,  find  the  intersection  point  (z£2, Xj2). 

•  Using  the  equation  for  the  plane  wave  ray  t?£  (z)  =  7nv(z  —  z£  2)  +  x^2  when 

mv  —  —cot (Gt  —  2o>£2)  —  cot (0o  —  2a>£2)  ,  0)i2  =  |tan_1(2C0Zi2)|  and  f(z),  find  the 
intersection  point  (z£1,  jr£1). 

•  Finally  use  equation  (36)  to  estimate  the  value  S£  for  this  point. 

Unfortunately  the  above  procedure  can  only  be  applied  when  the  location  of  the  grid  wire  is  appropriately 
inside  the  guide.  That  is  when  <f;-  >  <f2  in  Figure  15.  When  the  arc  length  is  smaller  than  <f2  we  have  to 
estimate  the  S)  in  a  different  way.  When  ^1<^£<f2we  need  to  use  a  modified  procedure,  and  finally 
when  <  fi,  yet  another  modified  process  has  to  be  employed. 


Figure  15.  Points  used  in  the  methodology 

3.8.  Computing  the  variable  angle  of  incidence  for  use  in  Marcuvitz’s  equation 

In  order  to  compute  the  grid  wire  diameters,  we  use  a  similar  methodology  as  before,  requiring 
the  equivalence  of  the  ideal  power  reflection  and  the  one  resulting  from  the  Marcuvitz  formula 
for  infinite  parallel  wire  grids  (see  our  1st  Quarterly  Report).  The  process  is  summarized  as 
follows. 

•  The  ideal  power  reflection  for  each  wire  location  on  the  arched  leaky  wall  is  now  computed  as 
follows  (now  using  the  computed  S£). 

Rp, ideal  Ofi)  =  exp(-5£a(^))  (37) 
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•  We  also  need  to  compute  the  power  reflection  using  the  generic  formula  Rpow  —  I F  [ 2  where 

F  =  ( Zeq  —  V)f(Zeq  +  1)  is  the  equivalent  circuit  reflection  coefficient.  The  formula  given  by 

Marcuvitz  (described  previously)  needs  to  be  iteratively  computed.  The  only  difference  in  the 
case  of  AAWSEA  compared  to  FAWSEA  is  that  the  angle  of  incidence  considered  in  Marcuvitz 
is  no  longer  constant,  varying  with  the  location  of  the  grid  wire  in  the  leaky  wall.  Figure  16  shows 
the  variation  of  this  angle  denoted  by  <plr  <p2, ...  in  the  figure.  When  the  angle  q>i  for  each  wire 
location  is  computed,  the  power  reflection  can  be  calculated  using  a  similar  Matlab  function. 

•  It  was  shown  in  section  3.4  (via  geometric  derivations)  that  q?i  =  —  60^.  Now,  employing  the 

perfect  compensation  assumption  =  60  for  i  =  1,  2, ...  that  will  become  q?i  =  80  —  0)^ 


The  required  angle  a)n  for  each  point  on  the  arched  wall  is  easily  computed  as  follows. 


.  \ 


ten  =  I  tan  V'Oii)!  =  tan  1 


Zil 


\JaZ  ~  Z‘J 

Now  substituting  zit  =  a  sin  in  equation  (38)  we  proceed  as  in 


(38) 


c On  =  tan  1 


Finally,  we  have 


(  a  sin  (^ 


=  tan-1  (tan  (~))  =  ~  (39) 


(40) 


Figure  16.  Variation  of  angle  of  incidence  on  the  leaky  wall. 

The  above  analyses  establish  most  of  the  framework  needed  for  an  (approximate)  AAWSEA  design 
theory,  to  be  tested  against  numerical  models.  We  are  currently  preparing  Matlab  scripts  and  functions 
that  incorporate  these  expressions,  and  which  will  allow  us  to  explicitly  compute  the  wire  sizes  and 
spacings  that  ought  to  go  with  the  subject  H-plane  curved  boundaries.  The  next  steps  will  be  to  compare 
the  predicted  wire  grill  configurations,  VSWR,  and  antenna  patterns  of  candidate  AAWSEA  configu¬ 
rations  to  our  more  established  FAWSEA  configurations,  initially  in  2D,  and  ultimately  in  3D  models  that 
include  dielectric  windows  and  multiple,  parallel,  finite-width  channels. 
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4.  DISCUSSION,  CONCLUSIONS,  AND  RECOMMENDATIONS 

Research  performed  during  this  9th  quarter  of  the  R&D  program  focused  upon  developing  a  theoretical 
framework/toolset  to  help  guide  the  AAWSEA  “design  recipe”  process,  proceeding  in  analogy  to  the 
theory  and  tools  we  established  previously  to  guide  the  design  of  FAWSEAs  and  CAWSEAs.  To  that 
end,  we  endeavored  to  expand  the  plane-wave  decomposition  perspective  into  a  more  general,  but  also 
more-approximate,  ray-tracing  technique.  More  work  is  needed  to  advance  and  validate  this  analysis, 
identify  its  limitations,  and  to  establish  its  effectiveness  in  accelerating  the  AAWSEA  design  process, 
which  (until  now)  has  been  based  mostly  on  geometrically  morphing  FAWSEA  designs  and  then  re¬ 
tuning  them  by  means  of  computationally-intensive  iterative  numerical  modeling. 

In  the  coming  quarter,  we  plan  to:  (1)  complete  the  investigation  (and  validation,  if  all  goes  well)  of  the 
generalized  ray  tracing  approach  (in  combination  with  the  previously-established  tools)  to  formulating 
AAWSEA  design  recipes;  and  (2)  make  significant  progress  toward  optimization  and  documentation  of 
“standard/  recommended”  designs  for  AAWSEA,  BAWSEA,  and  RAWSEA. 

As  always,  we  appreciate  ONR’s  continuing  support  for  this  R&D. 
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